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ABSTRACT 

The star formation rate - stellar mass relation (SFR-M*) and its evolution (i.e., the SFR main 
sequence) describes the growth rate of galaxies of a given stellar mass and at a given redshift. As- 
suming that present-day star- forming galaxies were always star-forming in the past, these growth 
rate observations can be integrated to calculate average star formation histories (SFHs). Using this 
Main Sequence Integration (MSI) approach, we trace present-day massive star-forming galaxies back 
to when they were 10 — 20% of their current stellar mass. The integration is robust throughout 
those epochs: the SFR data underpinning our calculations is consistent with the evolution of stellar 
mass density in this regime. Analytic approximations to these SFHs are provided. Integration-based 
results reaffirm previous suggestions that current star-forming galaxies formed virtually all of their 
stellar mass at z < 2. It follows that massive galaxies observed at z > 2 are not the typical progenitors 
of star-forming galaxies today. 

We also check MSI-based SFHs against those inferred from analysis of the fossil record - from 
spectral energy distributions (SEDs) of star-forming galaxies in the Sloan Digital Sky Survey, and color 
magnitude diagrams (CMDs) of resolved stars in dwarf irregular galaxies. Once stellar population age 
uncertainties are accounted for, the main sequence is in excellent agreement with SED-bascd SFHs 
(from VESPA). Extrapolating SFR main sequence observations to dwarf galaxies, we find differences 
between MSI results and SFHs from CMD analysis of ACS Nearby Galaxy Survey Treasury and Local 
Group galaxies. Resolved dwarfs appear to grow much slower than main sequence trends imply, and 
also slower than slightly higher mass SED-analyzed galaxies. This difference may signal problems 
with SFH determinations, but it may also signal a shift in star formation trends at the lowest stellar 
masses. 
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1. INTRODUCTION 

A history of star formation in star-forming galax- 
ies (SFGs) is central to the subject of galaxy evo- 
lution. To an extent, this history can be extracted 
from SFGs directly, by fitting stellar population mod- 
els to either resolved stars in color magnitude dia- 
grams (CMDs) o r to their spectral energy distribu- 
tions (SEDs; e.g.. iHeavens et al.ll2000t lAsari et alj|2007t 



Mateus et al.l 1200a lOcvirk et al 
20071: iCid Fernandes et al.l 120051 
Williams et al.l 120 111) 



" 120061: iToieiro et all 
Panter et al.l 120071: 

Unfortunately, an imperfect 
knowledge of stellar evolution, its slow pace in old pop- 
ulations and typical noise in observations, all conspire 
to produce broad uncertainty in ages that are estimated 
from this "fossil record" analysis. Average star formation 
histories (SFHs) can also be inferred from the tight star 
formation rate (SFR) - Af* relation because it describes a 
tight main sequence of star formation that connects SFG 
growth, stellar mass, and redshift. If present-day SFGs 
were always SFGs - always part of this main sequence in 
their past - then these growth rate observations can be 
integrated to calculate average SFHs. 

Data on the SFR main sequence is ripe for such 
exploration as multi-epoch, mass-binned SFR mea- 
surements, done through panchromatic surveys, are 



Oliver et al 


] [20ll iNoeske et al.l l2007al lElbaz et al.l 


2011L 2007 


: iMagdis et al.l 20101: iPannella et al. 20091: 


Ilbert et al. 


12010 


Karim et al.l 20111 Feulner et al.l 


20051 iBauer et al. 


120051 120111 IDrorv & Alvarez! 120081: 
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now report representative samples of galaxies that are re- 
sponsible for over 70% of star formation up to z w 2 (e.g., 
iKarim et ail 120111: IWuvts et alll2011bD . Moreover, the 
data suggests a tight SFR main sequence at all masses 
and redshifts probedZI- At z ~ 0.1, the SFR-M, relation 
is detected over more than four dec ades in stellar mass , 
from at least IO^Mq to 10 n M^ dWvder et al.l 120071: 
Brinchmann et all 120041: iSalim et alj 120071: lElbaz et al I 
20071 : iLee et al.l 1201 ID with approximately lognormal 
1(7 < 0.3 dex scatter. The small scatter in the rela- 
tion persists at z ~ 1 (INoeske et al.l l2007bl: lElbaz et al.l 
[2001 . z « 2 (iDaddi et al.l 120071: IRodighiero et al.ll20TTf r 
through z < 3 (|Elbaz et al.ll2011 ) and out to the earliest 
epochs at 4 < z < 7 ((Gonzalez et all 1201 ll where the 
relation may widen). 

If a SFG was always a SFG (i.e., formed stars through- 
out its history), then the SFR of that SFG's progenitors 
must appear in the scatter of the SFR main sequence 
at a given stellar mass and redshift. That SFR is the 
main component governing a SFG's mass growth; thus, 
starting with a present-day stellar mass, -M*o, a galaxy's 
growth rate from the main sequence (M*) can be inte- 
grated in lookback time to determine typical (z) and 
mass-dependent SFHs (see fJ3]for equations). Although, 
merging alters how stellar mass is distributed between 



Cowic & Bargcr (2008) should be noted as an outlier in this 



regard. 
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progenitors, we demonstrate in tJ3,l.ll that merging has 
a negligible impact on collective stellar mass evolution. 
If the assumption that SFGs were always SFGs is cor- 
rect, then the SFHs derived through this approach are 
entirely empirical. 

There is good reason to believe that SFGs formed 
stars throughout their evolution. At z ~ 0.1, a 
combination of characteristics that are not typically 
reversible can be used to reliably identify galaxies 
as quiescent. Quiescent galaxies h ave large stellar 
masses (e.g.. iBrinchmann et al.l 1 200-^. and dispersion- 
domi na ted kinematics (e.g~ [Roberts fc Havnesl 
1991 iKauffmann et all l200l IWuvts et al.1 1201 lat 



Schiminovich et al. and exist in r e gions of high 



galax y density (e.g.. iPeng et al.l 120101 : iSobral et al.l 
|2011[ ). At higher redshift, the quiescent fraction de- 
creases, meaning there are fewer quenched galaxies that 
could have been progenitors of pres ent-day SFGs (e.g., 
IBell et alJ[2007t IPannella et al.ll2009h . 

If the SFR main sequence reflected on-and-off bursts 
of star formation, then a quenched phase could be part 
of the duty cycle of SFGs; however, this possibility con- 
tradicts abundant evidence pointing to continuous star 
formation in observed SFGs. First, from a physical per- 
spective, the Sch midt-Kennicutt relation (Schmidt 1959; 
lKennicutt]|1998ft requires continuous consumption of gas 
in th e absence of violent disruption to gas disks. Fur- 
ther, INoeske et all (|2007b[) noted that the narrowness of 
the SFR main sequence itself suggests continuous evo- 
lution, and correlations between short and long term 
SFR indicators bear this out (e.g.. iQuintero et al.ll2~004t 
lYang et al.l [20081: iKriek et al.ll2011l) . The duty cycle in 
SFG progenitors is also constrained by the quiescent frac- 
tion, discussed above. Moreover, the distribution of star 
formation is bimodal, clearly se parating into two galaxy 
populations ([Wetzel et al.ll201ll ). Lastly, a low Sersic in- 
dex structural main sequence maps onto the SFR main 
sequence (see. IWuvts et al.l 12011a! iSchiminovich et ahl 
l2007t IKauffmann et al.l 12003ft . Thus, for bursts of star 
formation to be responsible for the SFR main sequence, 
they would also need to make galaxy surface brightness 
distributions less compact. 

Notably, many of these observations hold for the lowest 
stellar mass galaxies. In dwarfs, the dete ction of HI is 
always associated wit h star formation fe.g. JMeurer et ahl 
I2006t iLee et al.ll2011h . the dwarf SFR-M* relation is tight 
(|Lee et al.l 1201 ll ) . dwarfs also show a strong correlation 
betw een star formation over short and long timescales 
fe-g JLee et al.ll2009al: iMeurer et all 120061: iHunter et al.l 
120111) . and less t han 10% of isolated (central) dwarfs 
are quenched (see IWang et aJ J 120091: iHaines et all 120071: 
Uames et all 120081 : |Peng et al.ll201lF T 

Still, the assumption that present-day SFGs were al- 
ways part of the SFR main sequence may break down 
where quenching and SFR observations are not well con- 
strained. For example, at early epoch it is possible that 
quenched galaxies rejuvenated star formation following 
gas rich mergers. This type of scenario is important to 
bear in mind, so we highlight regimes where this Main Se- 
quence Integration (MSI) approach is extrapolated and 
observations are not yet representative or well under- 
stood (see ^2. 31 for further discussion). 

A number of other studies have also connected observa- 
tions of star formation at different epochs to the history 



of star formation in galaxies. INoeske et al.l (|2007aD was 
the first to use the SFR main sequence to infer SFHs, 
finding that, at z < 1, galaxies assigned a formation red- 
shift and exponentially decaying S FHs could match thei r 
SFR main sequence observations (jNoeske et al.ll2007b| ). 
Their "staged-r" model presents one way to account for 
the data, but there is no reason to restrict SFGs with 
that specific parameterization. A number of groups have 
made the more basic assumption that galaxies traced av- 
erage SFR-M* observations, thereby producing success- 
ful descriptions of the shape of the mass function, down- 
sizing, merger rates, the buildup of the red sequence, and 
rising SFHs in galaxies at z > 3 (see Drorv fc Alvarezl 
' 20081 IPeng et alj|2010HBell et al.ll2007l:lPeng et al.ll2011t 
Renzin i 2009; Papovich et al.ll2011l) . iLeitner fc Kravtsovl 



(|2011|) evolved SFGs backward starting at z = and 
showed how to incorporate mass loss self-consistently 
into integration of the SFR main sequence. The pur- 
pose of this paper is to calculate the quantitative impli- 
cations of the MSI approach using recent observations 
that may have recorded the early stages of present-day 
SFG growth. 

We will also confront MSI results with fossil record 
analysis that determines SFHs for individual galaxies 
by examining their SED s and CMDs. Simila rly, both 
IWalcher et all ([20081 ) and IWuvts et al.l (|2011fJ ) used the 
SFR main sequence to calibrate and improve SED- 
based SFHs, but the focus here is on understanding 
the different implications of MSI and the fossil record. 
The SED-based results in this paper draw on analy- 
sis of hundreds of thousands of gal axy spectra from 
the Sloan Digital Sky Survey (SDSS lYork et al.1 l2000l: 
IStrauss et al.|[2002[ ) but results are limited by well-known 
degeneracies, in terpretation of noisy data, and pitfalls in 
modeling (e.g.. IQcvirk et al.l 120061: ICid Fernandesl 120071: 
ICid Fernandes et al.l 120051: IRichards et all 120091) Using 
the SFR main sequence as a basis for understanding these 
SEDs allows a consistency check on MSI-based SFHs and 
an illustration of the effects of age errors on SED-based 
SFHs. 

CMD-based analysis is uniquely suited to probe early 
mass growth in the smallest galaxies. In this regime, 
MSI is no longer observationally constrained because the 
observations are only starting to probe representative 
samples of star formation beyond t he local universe at 
< 10 9 M o (e.g.. IWuvts et aljfeoiiah . By extrapolating 
SFR trends, we can compare MSI to local low mass dwarf 
SFHs and see how these smallest galaxies in the local vol- 
ume fit into the cosmological t rends (e.g., archaeological 
downsizing, iCowie et al.lll996h in galaxy formation. 

This paper can be divided into two main parts. First, 
in £j2] and in [JH we discuss SFR data and the SFR main 
sequence integration method, as well as the implications 
of MSI for galaxy growth. Second, in $4] we check for 
consistency between MSI and analysis of the fossil record. 
Specifically, we first describe recent SFR observations 
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and simple fits to SFR data in multiple surveys 
so as to estimate uncertainties. SFR data limita- 
tions are noted ( §2.3[) . both from survey incompleteness 
and inconsistency with the growth of stellar mass den- 
sity in the universe. We then review how multi-epoch 
SFR measurements imply the growth of stellar mass in 
SFGs (§2[). Accurate analytic approximations to SFHs as 
a function of stellar mass are provided in Appendix [AJ 
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Readers interested in the importance of scatter or merg- 
ers c an see ij3.ll Early growth in SFGs is described in 
i )3.2[ where results are compared to other scmicmpirical 
approaches (Abundance Matching and Halo Occupation 
Distributions). 

In the second half of the paper (2), we compare fos- 
sil record-based SFHs with MSI-based SFHs. de- 
scribes the samples we use to span the entire mass range 
of SFGs - from the ACS Nearby Galaxy Survey Trea- 
su ry (ANGST ) and L ocal Gro up (LG ) dwarfs analyzed 
in lWeisz et all ()2011aL hereafter IWllT). to SFHs from the 
versatile spectral analysis (VESPA; IToieiro et al.l 12007ft 
databasifl, which co mpiles their analysis o f the SDSS 
data release 7 (DR7, lAbazaiian et al.lT2009f ) . We com- 
pare MSI to fossil record results in §4.31 quantifying the 
importance of deconvolving errors for interpreting SED- 
based SFHs (a detailed discussion of age resolution can 
be found in Appendix [B| . Finally, in £14.41 we highlight 
the u nusual growth history of dwarf galaxies from the 
IWTTI LG+ANGST sample, using SDSS SED-based anal- 
ysis and MSI to place them in a cosmological context. 
A discussion of results in the context of simulations and 
other observations can be found in fJSJ We summarize 
and conclude in [J6] 

Throughout the paper a flat ACDM cosmo logy with 
D, Q = 1 - f2 A = 0.258 and h = 0.72 is used. A IChabrierl 
(|2003| ) initial mass function (IMF) will be assumed and 
all data is converted accordingly. 

2. MEASURING THE SFR SEQUENCE 

This section describes the measurement and reliability 
of the SFR data that we will use to constrain the growth 
of stellar mass in SFGs. 

2.1. SFR Data 

SFRs in external galaxies are always inferred from trac- 
ers of massive stars. These stars are short-lived and 
dominate the energy output of young stellar popula- 
tions. Their most obvious tracers are at short wave- 
lengths, including ionizing UV radiation as well as line 
emission from elements in the ionized HII regions (ob- 
servable in e. g. Ha and [OIL 1 ) surrounding O and early - 
type B-stars (jKennicuttl 119981: iBrinchmann et al.l 12001 . 
However, most of the energy from y oung stars - over 
80% in SFGs forming > 2M Pl yr~^ (IBuat et al.l l2010l: 
iTakeuchi et~aTll2010t IBothwell et al.ll2011l) - is absorbed 
by dust and re-radiated in the far infrared (FIR) spec- 
tral window. In fact, UV luminosity from young stars is 
subdominant to total IR luminosity (£tir, from 8/Ltm 
to 1000/im) in the local universe (z <~ 0.1) and con- 
tributes only marginally to the total bolometric flux of 
young stars observed in SFGs at higher redshifts (see 
Wuvts et alll2011aHTakeuchi et al.ll2005l: IBothwell et al.l 
20111) . Long wavelength observations, which can be more 
straightforwardly calibrated to reflect £tir, are therefore 
reliable tracers of star formation to early epochs. 

The Multi-band Imagi ng Photometer (MIP S) on the 
Spitzer Space Telescope (jWerner et al.l 120041 ) . has pro- 
vided a cruci a l wind ow onto dust SEDs in the mid-IR. 
iNoeske et all (|2007bD used MIPS-24/xm to measure the 
SFR main sequence and its scatter, thereby constraining 

2 S ee |http : //www-wf au.roe ■ ac .uk/vespa/| and IToieiro et al.l 
(2009). 



staged-T mo d els (INoeske et al.l l20~07al) . More recently, 
lOliver et al.l (|2010|, hereafter IO10J) analyzed the SFR 
main sequence over the Spitzer Wi de-area InfraRed Ex- 
traga lactic Legacy Survey ( SWIRE lLonsdale et al.ll2003l 
using the MIPS 70/xm and 160/im bandpasses. 
The lower sensitivity and resolution in these bandpasses 
necessitated stacking galaxies by stellar mass, but the 
longer wavelength observations avoided calibration un- 
certainties stemming from PAH emission fea tures^! , 
which redshift to_24um at z > 0. 8 (iFadda et al.1120061 ). 

IKarim et al.l IpOlll hereafter iKll " carried out the 
largest probe of the SFR main sequence to date by stack- 
ing radio observations from the Very Large Array (VLA) 
in the COSMOS field at 1.4 GHz between 0.2 < z < 3. 
Such radio SFR observations are dust model indepen- 
dent. However, the radio derived SFRs are only as good 
as their calibration to TIR luminosity: the physical origin 
of the tight empirical correlation between radio and TIR 
luminosity (|Belj l2003h is not trivially connected to star 
formation in widely varying interstellar medium (ISM) 
conditions. Fortunately, z = calibrations seem to be 
accurate to z = 2, but errors in the radio-derived SFRs 
may be domin ated by calibration uncertainty. 

Given iKllf s depth and large statistics , we use that 
study for our fiducial dataset. Together, IQIOI an d IKll! 
span the results reported in the literature. IQIOI results 
are therefore used to estimate analysis uncertainties as 
described in the following section. 

2.2. SFR Main Sequence Fits 

We integrate fits to the SFR main sequence data to de- 
rive MSI-based SFHs . The fits are done separately to the 
IkTTI (radio) and the [OlS (IR) surveys, thereby provid- 
ing an estimate of systematic uncertainties in analysis or 
calibration that are otherwise difficult to estimate. They 
will also be used to extrapolate measurements to regions 
of stellar mass and redshift where SFRs are not well con- 
strained (see §2.3.11 for further discussion). 

The SFR main sequence is commonly fit to survey data 
with a power law form in redshift and mass, 

TP(M*,z) = A 11 ( I - n ^-) (1 + zT. (1) 

Parameters a and f3 are then fit to data in redshift 
slices. Figure [1] shows the measured sSFR(z) at fixed 
M, = 10 11 M Q (which is characterized by An (l + z) a ; 
top left) and (3(z) (bottom left) for SFGs fromED and 
IQIOI The right panels of Figure [1] show data points from 
many recent surveys, demonstrating that me asure ments 
typically occupy the region between ECO and lOlOl 

3 Although, Herschel FIR observations showed that dust prop- 
erties are uniform across the luminosity function and to 2 = 1.5, 
such that old (e.g., (Chary & Elbaz 2001)) calibrations from L24/im 
to Ltir are accurate to < 0.2 dex at 2 < 1.5, only significantly 
overestimating SFRs in compact starbursts t hat d ominate sam- 
ples at 2 > 1.5 (identified by Elbaz et al. 2 0111; including in 
lElbaz et aflfe OlO; Rodighicro ct al. 2 010bl ; INordon et al.H2"OT0l; see 
also IBarro et al]|2011l ) . 

4 Herschel results show that the ratio of TIR to 1.4 GHz radio 
emission is c onsistent with a c onstant value w ithin the errors to 
2 » 1.5 ( see Hvis on ct al. 20U3 Figure 5, with IJarvis et all 120101; 
IBellH2003f) . At z < 2, results a re consistent with th e 2 = lijelB 
||2003|) calibration to ^ 0.2 dex USargent et al.H2010bl lat IMao et all 
l20Tll; [Bourne et al.ll2oTll. Figure 13). 
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Table 1 

Fit Parameters for ip(M t ,z). 
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1.15 


Kll 
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Note. — Small formal error bars have been omitted, meaningful uncertainties 
are given by discrepancies between different survey results. 

*Fits to specific SFR and /3sfg are based on the redshift-binned data shown in 
Figure [TJ (not binned by mass). 



The (3{z) parameter in Figure [TJ (bottom) shows the 
extent to which galaxies of different masses assemble at 
different rates. A negative value means more massive 
galaxies have smaller specific SFR (sSFR = SFR/Af*) 
and, therefore, must have formed a larger fraction of their 
stellar mass at high redshift. This is the "archaeologi- 
cal" form of downsizing, and is the only form of down- 
sizing referred to throughout this paper. Some downsiz- 
ing is evident from sSFRs - a process also implied by 
mass function evolution ( e.g. iPerez-Gonzalez et alJl2008l : 
iConrov fc Wec hsler 2009 ). Th e value of (3 and its evolu- 
tion vary by survey, but Kill showed that this variation 
can be largely explained by fluctuations at the edge of 
mass ranges, and aggressive color selection within the 
SFR main sequencCl Despite the use of a wide vari- 
ety of different methods and calibrations, there is rough 
agreement in sSFR observations, in the sense that trends 
are larger than the scatter between different studied 

In practice, to roughly characterize observations of the 
SFR main sequence report ed by various studies, we fit 
three different forms to the lKUl and lOld SFR data: 

1. A power law form ("pow") that follows Eq. [I] with 
constant a and /3. The IK 111 fit is our fiducial 
model. The iQICil fit includes a break to a constant 
ip{M*) at z > 2. 



A power law times exponential form ("powexp") 
that follows Eq. [2] with constant a, (3, and 77, 



4(M„z) = A 1 



M 



10 n A/« 



/3+1 



(l + z) a e' )z . (2) 



We thereby account for curvature in Figure [JJ 
(left), and the form allows a formally consistent 
fit to the [Kll] SFG data (reduced-* 2 = 1.15 and 
n dof = 44). 

3. A non-parameteric form that linearly interpolates 
(in 1 + 2) between measurements at different 

5 Meaningful selection of "star-forming" galaxies is not ambigu- 
ous. At all masses the galaxy bimodality ide ntifies quenched and 
star-forming ga laxies based on str ucture (see, IWuvts ct al. l2011aT ) 
and SFR (e.g., [Wetzel ct al. 201 Jl) . If the selection criteria accu- 
rately reflect the galaxy bimodality (see discussion in IK11T1 , the 
specific details should not affect the recovered median SFR main 
sequence. 

6 Note error bars are generally underestimated (e.g., left panels 
of Figure [TJ. 



redshifts ("interp"). Outside of the data range, 
cons tant and linear extrapolations are used for 
IQlOl and IK 111 respectively. 

Parameters of the po wer la w and exponential fits are 
reported in Table [TJ For IK 111 data, we derived the max- 
imum likelihood a, j3 and 77 from their average stacked 
SFR results and errors, binne d in both mass and redshift, 
from their Table 3. For lOlOt the fits are done separately 
for the sSFR at M* = 10 n M Q and /?, based on the points 
and errors plotted in Figure [TJ Fits are plotted as lines 
in the left panels of Figure [TJ 

2.3. Reliability of Results 

Integration of SFRs to derive SFHs is limited by the 
availability of SFR main sequence data and the reliability 
of that data. 

2.3.1. Completeness Limits and Extrapolations 

The stellar mass at which SFGs are well represented 
in a dataset is a "com pleteness" threshold. Our fidu- 
cial dataset from IK 111 reports their 95% statistical com- 
pleteness threshold - the stellar mass at which the VLA- 
COSMOS flux limits and source confusion still allow 
> 95% of SFGs to be stacked in order to measure SFRs 
at a g iven redshift. Strictly speaking completeness lim- 
its in IK 111 make it impossible to derive star formation 
histories for dwarfs, and limits firm results to massive 
(> 1O 1O M ) disks. 

However, these limits are conservative, and do not 
mean that extrapolations arc unreliable. On the con- 
trary, there are no clear indications that the tight SFR- 
M* relation falls apart at any mass or redshift. Locally, 
SFR main sequenc e trends can be traced a c ross the en- 
tire mass function (iBrinchmann et al J 120041 : iSalim et al.1 
120071 : iWvder et al J 12007( 1 . with sSFRs rising t oward low 
masses and staying high in dwarf galaxies (e.g. lLee et al.l 
I2009U IJanTes et alj|2008t iLee et alJlMTJ) . Furthermore, 
the s equence is also in pla ce at low masses at high red- 
shift (jGonzalez et al.ll201lH . More recently. IWuvts et al.l 
(|2011all mapped the sequence to 10 8 M© and out to 
z = 2.5, finding no obvious change in galaxy properties 
at any mass or redshift. 

Still, data remains limited and hints of a flattening of 
sSFRs in low mass galaxies and at high redshift have been 
reported (see IK 111 for further discussion). Regions where 
MSI-based star formation histories rely on extrapolated 
trends are therefore noted with diagonal stripes in figures 
below. 
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Figure 1. Top panels: the sSFR of 10 Mq galaxies as a function o f reds hift. Bottom panels : fits to the powerlaw slope of the 
SFR-Af* relation as a function of rcdshifts. Left panels: lines show fits to lKlll (magenta squares) and lOldl (cyan triangle) SFR data (see 
the legend and numbered text for forms). Right panels: the shaded region marks the area enclosed by our fits (lines) from the left panels. 
The vertical dotted lines mark z eTT from Figu re |21 Points are a compilation of fi ts to SFR main sequence in re dshift b ins from a number of 
groups (hollow triangles: Nocskc et al. 2007b; hollow squares: Elbaz ct al. 2007; pentagons: Salim ct al. 2007; circles: Pannclla ct al. 2009; 
crosses: [Daddi ct al. 2007; three-line cross: IDun nc ct al. 2009; asterisks: [Rodighicro ct al. 2010bl). 



6 



2.3.2. SFRs at z > 1 and Consistency with Stellar Mass 

Density 

A number of studies (e . g.. iHopkins fe Beacoml 1200 
Conroy fc Wechslerl |200a iWilkins et all l2008bt iDav 



0.14 



z = 0.1 



0.5 



2 3 5 



20081 iRudnick et all 12009 ) have reported that the in- 
tegral of the cosmic SFR density from z > 2 ex- 
ceeds the growth in stellar mass density relative to 
what is observed (modulo mass loss). Since measure- 
ments of stellar mass density are not subject to the 
same extrapolation from high mass stars to the bulk 
of stellar mass , the discrepancy co uld i mply an evolv- 
ing IM F (e.g. iHopkins fc Beacoml 120061: IWilkins et all 
2008albT: iFardal et all I2007t iPerez-Gonzalez et all 120081: 
van Dokkuml 120081 : iDavel 120081 : iGunawardhana et all 



20111: IBorch et alJ 12006ft ." or a broader failure to prop- 



erly model the SEDs of these denser, mor e active, lower 
meta l licity, high-redshift galaxies (e.g., iNordon et all 
120101: lArnouts et all 120071 : altho ugh these are increas- 
ingly constrained by Herschel, see (Elbaz et alJl20Tl . 

In Figure[3J we compare the cosmic SFH from Kill with 
predictions for the cosmic SFH b ased on the buildup o f 
stellar mass from a compilatio n bvlWilkins et all (|2008aD . 
Adjusting for their IMF, the IK 111 points are consistent 
with the 1-cr errors on stellar mass growth predictions 
below z e rr = 1.6 (marked by the vertical dashed line). 
SFR measurements at > z orr are discrepant by > 3o\ 
The stellar populations that are measured, and the cor- 
responding MSI results, cannot be reliably interpreted 
at z > z Grr . MSI-based findings at z > z orr arc therefore 
noted with dotted shading in figures below. 

In the context of past cosmic SFR measurements. | K11 



findin gs are on the low end of th e IHopkins fc Beacom 



(|2006l ) SFR compilation (see ECD Figure 11), and so in- 
dicate that the IK 111 S FR main sequence is an outlier. On 
the other hand. lKllI is consistent with both past IR data 
and the growth of stellar mass density in the universe at 
z < z rrT . Moreover, the total star formation observed by 
IKIII since z err i s similar to the total median star forma- 
tion inferred bv lWilkins et all (|2008aD from the evolution 
of stellar mass density; applying a corr ection for a hypo- 
thetical evolving IMF would drive lKlll measurements out 
of agreem ent w ith SFR data. 

While lOlOl do not measure a cosmic SFH, 
iRodighiero et all (|2010al ) included FIR measure- 
ments of the same SWI RE fi eld in their measurements 
of the cosmic SFH, and Kill report that these measure- 
ments are consist ent w ith their findings (within large 
error bars). The lOlOl measurements can therefore be 
considered valid in a similar redshift range, but with 
more stringent completeness limits because of their 
shallower SFR observations. 

3. THE MSI APPROACH: FROM STAR FORMATION 
RATES TO STAR FORMATION HISTORIES 

With a reliable SFR main sequence in hand at z < 1.6, 
observations can be synthesized through MSI to calculate 
SFHs, <E>(i). The evolution of stellar mass in a given 
galaxy is described by a growth term from the observed 
SFR and a loss term from the galaxy-wide stellar mass 
loss rate, 

M*(i) = i)(M*,z)-n(t). (3) 

Galaxy-wide mass loss is given by an integral over the 
fractional mass lost rate (/ m i) from each single age stellar 




4 6 8 
redshift 



Figure 2. Cosmic SFR from the IK11I SFR sequence data com- 
pared to the cosmic SFR ex p ected g iven the growth of stellar mass 
measured by Wilkins ct al. (20Q8a|). The gray region marks the 
IWilkin s ct al. (2008a) la and 3 <7 uncertainty region adjusted from 
their IMF to a lchabrierl H2003 1 IMF by adding 0.045 dex based on 
the difference in b olometric flux for a bursting pop ulation (calcu- 
lated in PEGASE.2 IFioc fc Rocca-Volm crangc 1993) ■ The vertical 
line marks the redshift above which SFRs are not well understood 

(Zcrr = 1-6). 



population (SSP) - calculated u sing the flexible stellar 
populations synthe sis code (v2.0 iConrov fc Gunnl 120101 : 
IConrov et aD lOQS or see iLeitner fc Kravtsovll2011l for 

fits), 



»(t) 



*(0/ml(* - t')dt'. 



(4) 



to 



The SFH of that galaxy is then, 

$(i)=V(M*(t),«), (5) 

where ip is a function of the evolving mass. $(£) and 
M*(£) can be thought of as the lagrangian coordinate 
tracking a given galaxy, while ip, the SFR main sequence, 
is the SFR at a fixed Eulerian coordinate in the space of 
stellar mass and redshift. 

For our fiducial choice, galaxies start with boundary 
condition M*o = M*(z bs), and trace t he me dian of the 
SFR main sequence given by Eq. [Tj and IK 111 parameters 
(in Table [TJ as they move back in time and down in 
stellar mass. Since $(£) is needed to calculate mass loss, 
but mass loss is needed to solve Eq [3l a self-consistent 
solutio n for M* (t) requires iteration. ILeitner fc Kravtsovl 
(|2011D described the simple iteration procedure that is 
used here to converge on a self-consistent $(f) and 3?(£). 

Figure [3] shows mass growth for galaxies including the 
spread from variation between different SFR main se- 
quence observations (described by fits in £|2.2|) . Ap- 
pendix reports analytic approximations. MSI-based 
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SFHs decay almost exponentially after early buildup, and 
are thus i n general agreement with staged- t models as- 
sumed bv lNoeske et all (|2007al) . 

3.1. Assessing MSI Complications 

In this section we address two issues - mergers and 
scatter - that were ignored by linking galaxies across 
SFR observations solely using the median SFR and re- 
sulting stellar mass evolution. We find that these issues 
are probably less important than variations between ob- 
servations of the SFR main sequence, so readers inter- 
ested in results can skip to £13.21 

3.1.1. Mergers 

MSI is almost insensitive to mergers, both because sS- 
FRs of SFGs are not strongly dependent on mass and 
because mergers are not common in SFGs at z < 2. 
Starting from a M*(z b s ), galaxies disassemble accord- 
ing to their growth rate, M*, with increasing lookback 
time. If at some z sp ii t there is a merger of ratio /x mra t, 
the galaxy gets split into two progenitors, both of which 
continue to contribute to the SFH of the final galaxy. 
Since sSFRs are mass dependent, the same galaxy grows 
at a different rate if it is split into progenitors. The ratio 
of growth in a split to unsplit galaxy is 



7SFR 



1 + Mmrat 



(1 + Mmrat) 



1+/3 ■ 



(6) 



Observations indicate that galaxies downsize (/? < 0, see 
fc|2.2[) . so low mass galaxies grow faster and 7sfr > 1. 
Deviations are maximized for rare equal mass mergers, 
and 7s fr < 1-3. 

To confirm that merging is not important to stellar 
mass growth, we couple MSI with merger rate estimates, 
thereby turning MSI into a empirical generator of merger 
trees. The merger rate per galaxy, v(M*(t), z, /U mra t)j is 
tabulated from the semi empirical Halo Occup ation Dis- 
tribution based results of lHopkins et al.l ([20101 ) with code 
that they provided- The (dis-) assembly is then realized 
with 1,000 Monte Carlo simulations, where each galaxy 
and its progenitors evolve according to both their ob- 
served SFR, t/j(M*,t), as well as by splitting their stellar 
mass into progenitor systems at each timestep with prob- 
ability i/(M*(i), z, ^ mra t)<^EI SFRs and stellar masses are 
averaged at each timestep and the resulting SFH of the 
average galaxy due to merging is plotted in Figure 01 
along with the same fiducial simulation without merg- 
ers. Merging results in very slightly younger SFGs. 

Note that this is in no way a physical statement about 
merger rates or satellites. It is an empirical observation 
that splitting a SFG into its typical constituent progen- 
itors should not (on average) substantially change the 
collective growth history of the galaxy. 

3.1.2. Scatter in SFR-M* and Galaxy Environment 

In reality, galaxy SFRs scatter around the me- 
dian SFR-M* value, with approximately lognormal 
0.3 dex scatter. (lElbaz et al J 120071: iNoeske et aLll2007U 

7 http: //www. cf a. harvard . edu/-phopkins/Site/mergercalc .html 

8 This is similar to the method used in Hopkins et al. (2009) to 
explore bulge formation. 



lElbaz et "all 120111 ). Taking Awsfr) to be the devia- 
tion of a given galaxy from the median relation in some 
timestep, there are two extreme scenarios for populat- 
ing that scatter: the Ai og (gFR) can be constant across 
timcstcps, or Ai og (gFR) can be totally uncorrelated be- 
tween timesteps. In the correlation coefficient c f = 
case (no covariance), average galaxies form stars at the 
average SFR inside of the lognormal scatter (about the 
median, ip). Thus the Ct=0 case closely resembles the 
fiducial case (with the mean SFR replacing the median). 
In the case where Ct = 1, active SFGs with Ai og (gFR) > 
are always active and hence formed a shorter time ago. 
Correspondingly, less active SFGs (A log ( S FR) < 0) are 
always less active and therefore formed at earlier epochs. 
When averaged, the result is a smeared out SFH as 
demonstrated in Figure [4] 

Observations that parse the breadth of the SFR-M* 
relation for SFGs could implicate one of these scenar- 
ios. Insofar as galaxy environment affects accretion his- 
tory, a second-order correlation of Ai og ( SF R) with galaxy 
envi ronment would be t he most likely cause of c t ^ 1 
(see iDutton et al.1 120101 ). A number of groups have 
found persistent c orrelations between environment and 
SFG fraction (e.g . . iKauffmann et a~ 120041: ICooper et al 

ll 



20101 iPatel et al.l [20111: iSobral eta! 



lHaines et al 



2007t) . However, focusing exclusively on SFGs, there are 



no conclusive indications that A log (g F R) is signifi cantly 
driven by environment. At z ~ 1, ISobral et al.l ([201 1[ ) 
found that SFGs form stars ~ 50% faster in group en- 
vironment^ relative to more isolated galaxies of M, < 

4 x 1O 1O M (consistent with lElbaz et al.1 12007ft Mean - 
while, at z sa 0.1, findings are mixed: iPeng et al.l (|2010f ) 
discerned no difference b etween SFRs at fixe d mass and 
varied environment, but lHaines et al.l (j2007|) claimed a 
significant detection of an inverse Ai og (gpR)-density re- 
lation in low mass SFGs and a margina l inverse rela- 
tion spanning ~ 0.3 dex is also apparent in lPopesso et al.l 
([201 ll see Figure 15) from their low O.lMpc -2 < £ < 

5 Mpc~ 2 to intermediate 5 Mpc" 2 < £ < 30 Mpc~ 2 den- 
sity environments. Clearly it is plausible that c t = 0, but 
there is also not enough information to conclude c t ^ 
- that either uniform trends in SFG fraction don't carry 
over to the SFR- A/* relation (c t = 1), or that SFGs en- 
counter an environmental inversion. 

To check for the impact of a (hypothetical) environ- 
mental inversion, we re-simulate average MSI SFHs tak- 
ing the Ct = 1 model from the previous section, and 
impose an inversion in the amplitude of the scatter for 
each Monte Carlo simulation at the Zinvcrt = 1 time-step. 
Figure S] shows the average of the realizations. This toy 
model introduces a sharp feature in the average SFH be- 
cause z < 1 active isolated galaxies have burned them- 
selves out, leaving a population dominated by relatively 
passive systems that are relatively active at z > invert- 
The resulting average mass growth and SFHs are similar 
to the fiducial model. 

Other secondary correlations with Ai og ( SFR ) may fur- 

9 Merging SFGs in clusters form stars rapidly, but these are rare, 
and since such bursts occur in regions with low SFG fractions to 
begin with, they may quench star forma tion directly (a con clu- 
sion also inferred from st ructural data by Wuyts ct al. 20 11al and 
ISchimi novich ct al. 200?]). As p art of an ext ended tail in the oth- 
erwise lognormal sequence Csee. lElbaz "c t al. 201]]), they arc likely 
not relevant for the history of present-day SFGs. 
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lookback time [Gyr] lookback time [Gyr] 

Figure 3. The SFHs (left) and mass growth (right) from MSI, in galaxies of M t0 = 10 s Af Q , 1O 9 M , 1O 1O M and W^Mq. Shaded 
regio ns show the variation between SFR main sequence observations and thick lines are fiducial power law SFR main sequence fit to 
IK11I data. Results are sho wn whe re SFR main sequence data is robust and complete (cross hatching), extrapolated (diagonal stripes; see 
completeness discussion in i|2,3.1H and compromised by observational uncertainties at z > z crr (dotted). 



ther constrain the SFHs of SFGs. Properties, such as 
galaxy clus t ering , circular velo city, m orphology (e.g., 
lElbaz et al.l 12007): IWuvts et all I2011al ). the evolution 
SFG number counts, and merging or gas physics, may 
be used to identify the position of galaxies in the scatter 
of the SFR main sequence. Each place important pri- 
ors on SFG growth (see, e.g.. iSchiminovich et al.l 120071 : 
iBoissier et atll2010t iBouche et al.ll2010f ). However, the 
focus here will be on the implications of consensus SFR 
observations in a straightforward empirical sense, rather 
than on embedding the integration in a more physics- 
driven semi-analytic model. These observations are 
enough, by themselves, to generate interesting and easily 
interpreted constraints on galaxy growth. 

3.2. The Delayed Assembly of Star-Forming Galaxies 

Figure [3] indicates that measuring the early stages of 
present-day SFG evolution does not require observations 
of the dawn of galaxy formation. Galaxies observed in 
the SFR main sequence at high redshifts (e.g., z > 3) are 
mostly quiescent and massive by z = 0, while present-day 
SFGs appear to have grown most of their mass starting 
at relatively low redshift (1 < z < 2). 

3.2.1. Delayed Growth and Relevance to Bulge Formation 

SFR main sequence measurements from IK 111 are deep 
enough that MSI can follow the growth of stellar mass 
back to a time when typical SFGs were almost a tenth of 
their current size. Figure [5] (left) shows the redshift after 
which MSI indicates average SFGs formed 15% of their 
stars (z 15 %). The thin lines in the figure show the ef- 
fect of a factor of two increase or decrease in the fiducial 
SFR main sequence; this can be interpreted as illustrat- 



ing large systematic uncertainties, or a maximal estimate 
of the 68% scatter in SFHs (i.e., in the Ct = 1 case from 
§3.1.21 although ignoring z err ). 

We focus on z 15 % specifically because, for 10 10 M Q < 
M*o < 5x 10 10 , constraints on z 15 % are still data-driven; 
z 15 % results are affected by neither incompleteness, nor 
disagreement with mass function determinations (z orr 
from Figure [2]). This regime also bears particular rel- 
evance for disk galaxy structure and bulge formation in 
simulations since median observed bulge mass fractions 
are also 15% (see N5.1l for further discussion). Moreover, 
focusing on this mass fraction circumvents complicating 
factors related to bulge formation (e.g., rejuvenated star 
formation from mergers). 

The figure indicates that the peak of star formation 
in SFGs is delayed relative to high-rcdshift (z > 2) star 
formers: 85% of stellar mass in many SFGs appears to 
form after z « 1.5. 

3.2.2. Model Comparisons at z = 1 

Complementary scmicmpirical methods, based on the 
galaxy-dark matter connection and simulations, also 
expect delayed formation of SFGs. Figure [5] (right) 
shows the fraction of mass formed in galaxies by 
z = 1, including result s from Abundance Matching 
(|Conrov k, Wechslerll2009l data from their Fig ure 3) and 
Halo Occupation analysis (jZheng et al.ll2007bT data from 
their Figure 7). These methods paint galaxies onto dark 
matter halos in simulations at different epochs accord- 
ing to mass (Abundance Matching) or clustering (Halo 
Occupation), and then link galaxies across epochs us- 
ing average dark matter accretion rates onto those dark 
matter halos during the intervening time. Note that 
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Figure 4. Changes to MSI-based SFHs (left) and stellar mass growth (right) due to mergers and scatter in SFR-Af* for SFGs of 
Af»o = 10 10 ' 5 Mq. The gray solid region encompasses the spread between MSI based on IK11I and IQlOl data. The blue line shows the 
fiducial model to which other lines should be compared. The green dotted lines show averaged MSI merger simulations. The dashed lines 
show the average of Monte Carlo realizations of the lognormal SFR PDF: black short dashed shows a constant Ai og (sFR) across time in 
a given galaxy and the red long-dashed line shows a constant Ai og (gpjy that inverts at 2 = 1. The vertical dotted line is at z CTT from 
Figure [2] 



these methods arc applied to the galaxy population as 
a whole, rather than being restricted to SFGs; however, 
at M„ ~2x 10 10 M Q galaxies transition to being mostly 
star-forming so the comparison is useful. 

Both Abundance Matching and MSI predict lower 
mass galaxies formed a large portion of their mass af- 
ter z ps 1, while galaxies of ~ 10 u Mq form most 
of th eir mass at higher r edshif ts. Unexpectedly, how- 
ever, IConrov fc Wechslerl (120091) find that even less stel- 
lar mass was in place in z = 1 progenitors than MSI in- 
fers - the addition of quiescent galaxies sh ould increase 
the fra ction of stars that form early. The IZheng et al.l 
(|2007bf ) Halo Occupation analysis predicts that even 
smaller stellar masses were formed before z = 1 rela- 
tive to Abundance Matching. For now, systematic er- 
rors in the halo-stellar mass relation can account for 
the 0.3 dcx tension between different predictions (e.g., 
from errors on stella r masses and sample va riance, see 
iBehroozi et~"ai1 120101: iLeauthaud et all 120111 ). Indeed, 
Figure 12 of IBehroozi et al.l (|2010l ) explicitly shows that 
IZheng et al.l (|2007b[ ) halos host fewer stars than most 
other determinations predict at z = 1. With tighter con- 
straints, differences from halo-based predictions might be 
used to infer star formation duty cycles or the possibility 
of a quenched phase. 

4. COMPARISON OF MSI WITH FOSSIL RECORD 
ANALYSIS 



SFHs are more traditionally derived through analysis 
of a galaxy's fossil record. This analysis uses either a 
galaxy's SED or, in the few local cases where individ- 
ual stars can be resolved, stellar CMDs. Given an SED 
F\, a stellar census &(t agc ) can be calculated by invert- 
ing F x = J $(t agc )S\(t agc , Z)dt agc , where S\{t,Z) is the 
SED produced by a SSP in a population synthesis model 
that includes the metallicity (Z) dependent stellar evo- 
lution tracks, an IMF, and dust extinction. CMD-based 
SFHs are similarly derived, but in this case stellar tracks 
are fit directly instead of being summed over. 

Unfortunately, this inverse problem is strongly ill- 
conditioned; low levels of noise in Fx can result 
in widely varying ages for stellar population, with 
serious i mplications for interpretation of SED-based 
SFHs CseelMoultaka et alll2004t iMoultaka fc Peladl200tt 
lOcvirk et al.l 120061; iCid Fernandes et al.l 120051 120011) . 
Moreover, large statistics can insidiously shrink error 
bars without improving recoverable resolution. Still, fos- 
sil record analysis provides important clues to the growth 
of mass in galaxies, and is singularly important for mea- 
suring the growth of dwarf galaxies where SFR main se- 
quence data is sparse. 

4.1. Fossil Record Data: Samples from SDSS and 
LG+ ANGST Dwarfs 

In this section, we briefly describe the fossil record sam- 
ples and analysis methods used. Table [2] summarizes the 
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Figure 5. Left panel : the reds hift by which 15% of stellar mass formed in all progenitors of a galaxy from MSI. The shading encompasses 
variation between Kll and O10 SFR main sequence observations. Striping indicates the reliability of the results as desc ribed in Figure l3l 
The dashed lines show, for reference, the effect of large scatter (or error) in 215% calculated by adding ±0.3 dcx to the Kll (thick blue) 
SFR main sequence. Right panel: the average stellar mass fraction formed by z = 1 in all progenitors as a function of z = stellar mass 
from MSI (blue as above), Abundance Matching (red dashed), and Halo Occupation modeling (green dotted). 



Table 2 

Galaxy Samples 



Mass Range [M ] 


N gal 


{M,(* obs )[Af Q ]> 


(Zobs) 


<tlb[Gyr]> 


Source 


M* <> 10 8 


32 


4.0 x 10 7 


D < 4Mpc 


~ 


LG+ANGST Wll dlrr (CMDs) 


10 s - io 8 - 5 


1,115 


2.1 x 10 8 


0.014 


0.19 


SDSSspg (SEDs) 


10 9 - io 9 - 5 


12, 156 


2.1 x 10 9 


0.032 


0.43 


SDSSspg (SEDs) 


10 9 - 10 10 


39, 390 


5.0 x 10 9 


0.042 


0.57 


SDSSspg (SEDs) 


10 io _ 10 io.S 


50, 071 


2.0 x 10 10 


0.070 


0.92 


SDSSsfg (SEDs) 


10 io.5 _ 10 n 


55, 780 


5.8 x 10 10 


0.098 


1.26 


SDSSsfg (SEDs) 


10" - io 11 - 5 


29, 187 


1.7 x 10 11 


0.12 


1.61 


SDSSsfg (SEDs) 



relevant properties of the sub-samples from which SFHs 
are derived. 

SED-based SFHs (for galaxies with M * > 10 8 M e ) 
are drawn from the VESPA database^ (jToieiro et al.l 
120071 120091 ) . which stores SFHs from VESPA applied 
to the spectroscopic sample of the of the SDSS DR7. 
SFGs are selected from the main galaxy sample (MGS 
iStrauss et all |2002[) using emission line measurements 
from the general purpose SDSS pipeline, in order to 
exclude both quiescent galaxies and active galactic nu- 
clei(AGNs). SFGs are required to have Ha emission de- 
tected at > 3er sig nificance. If Ha, H/3, [Nil] (A6584A), 
[OIII] (A5007A) are all detecte d at > 3cr, then we ex clude 
AGN with a cut in the BPT (| Bald win et al.lll981l ) dia- 
gram from|KM^anEitin (HM Eq. 1), log (^gjf) < 



0.61 



\ogfm 



1 http: //www-wf au. roe . ac.uk/vespa/ 



^^j^r 0.05 J + 1.3, and the requirement 

that NII/Ha < 0.6; these cuts are in the spirit of 

IBrinchmann et aT] (|2004l) . No efforts have been made 
to account for galaxy selection, but we note that harsher 
color cuts (e.g., U — R < 1.5) selecting more active SFGs 
had no effect on our results. Placing a low S/N thresh- 
old (median S/N> 10 in all bands with unmasked VESPA 
(S/N) > 20) also made little difference to SFHs. 

The total SFG sample includes about 175,000 galax- 
ies that are split predominantly into three VESPA mea- 
sured stellar mass bins of lO 9 - lo M ,lO lo - ia5 Af o and 
10 10 - 5 - n M o . Two smaller bins, at very low (10 8 " 8 - 5 M Q ) 
and high (1 11 ~ 11 5 M ( t ) ) mass, are also recorded. 
iToieiro et at] (|2007l) have shown that VESPA SFHs 
are consistent with SEP analysis from textttMOPED 
([Heavens et al.l 120001 ) , but VESPA also somewhat miti- 
gates over-fitting pitfalls by limiting parameters recov- 
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ered as advocated by lOcvirk et al.l (I2006T ) . The sample 
therefore represents the best SFH analysis of an unsur- 
passed number of galaxy fossil record observations. 

For lo wer ma ss dwarf galaxies, we use the SFH compi- 
lation bv lwTTlfrom the ANGST (iDalcanton et al.l l200l 
measured by IWeisz et al.l (|2011bD , and from the LG 
' Mated 119981 Table 1). with CMD s from iDolphin et all 



20051 ) and iHoltzman et~aTl (|2006t ). The sample is vol- 
ume limited to D < 4Mpc. Individual stars were re- 
solved and fit with stellar tracks in a CMD. Present-day 
type and color has little bea ring on past SFH beyond 
the last 1-2 Gyr (jWeisz et alJl2011bl iaf). but we will con- 
fine our discussion to the relati vely is olated LG+ANGST 
dwarf irregular (dlrr) sample (|Wlll ) as a low- ma ss ex- 
tension of our normal star-for ming sam p le. The lWTTl 
analysis was performed with a iSalpeterl ()1955[ ) IMF, so 
s tellar mas s es are multiplied by — 0.25dex to convert to 
a lChabrierl ()2003[ ) IMF. Their detection limits make IMF 
difference s unim portant for cumulative star formation 
plots fsee lwTTf ). 



4.2. Uncertainties in Fossil Record Analysis 

To account for age errors, we note that SSPs change 
logarithmically with age intervals and therefore convolve 
MSI results with a Gaussian filter in log-age. In or- 
der to choose a filter width, we note that the SFR 
weighted average bin size in VESPA is ~ 0.4 dex and co- 
variance between adjacent bins fur ther reduces resolu- 
tion. For much higher S/N SEDs, lOcvirk et ail (j2006[ ) 
found that disentangling confounding degeneracies, cou- 
pled with meager SSP differences, lead to age resolution 
(full width at half-maximum , A ag0 ) limits of > 0.8 dcx. 
Below wc present results smoothed with A agc = 1.0 dex, 
but findings are unchanged at > 10 9 Mq assuming any 
0.5 dex < A ago < 1.0 dex (see Appendix IB1 for further 
d iscuss ion) . 

I Will analysis reports broad age bins, such that covari- 
ance in simulated data result in uncertainties that are 
within measurement errors; but, again, this covariance 
cannot be miti gated by statistics (see their Appendix A 
and Figure 13. IWeisz et al"1l2011br) . Age errors in CMD 
analy sis can be smaller than in SEDs (see iGallart et al.l 
120051 for a discussion), but our CMD-related results be- 
low are qualitatively unchanged by the error model, so 
we will use the same smoothing filter for all comparisons. 

In addition to resolution, there are substantial system- 
atic uncertainties in stellar populati on synthesis (SPS) 
models. VESPA analys is, usin g both | Bruzual k, Charlotl 
(here after BC03; l20q3l ) and the iMarastonl (hereafter M05; 
120051 ) SPS models, is shown in plots. Those differences 
render Poisson error on the mean (of order the point 
size ), irre levant for all of the SDSS samples in Table [2l 
For IWI il dlrrs, statistical uncertainty is similar to un- 
certainty between stellar evolution tracks (see their Fig- 
ure 3). In both cases, error estimates are only qualitative 
in the sense that they illustrate the difference between 
models and not the uncertainty in the model parame- 
ters. 



Figure p shows the average cumulative star formation 
fractioiCJfor VESPA analyzed SDSS SFGs and MSI, di- 
vided into panels for different stellar mass bins. MSI 
models were initialized at the average mass and redshift 
of the SEP comparison sa mple (Table Et. Black points 
mark IBC03I (squares) and IM05I (circles) SPS analysis. 
These models show substantial differences that are also 
visible in the (logarithmic) plots of pre vious studies (e.g., 
iPanter et aLll2007tlToieiro et al.ll2009D . likely originating 
in their different treatment of thermally-pulsing asymp- 
totic giant branch (TP-AGB) stars. The se discrepancies 



remain the subject of active debate (e.g-JMarasto 



Conrov fc Gunnll20Tol: IMarastonll2011l: IKriek et al.l 
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4.3. Consistency with SDSS SEDs 



MacArthur et alj|201 

Focusing on the VESPA results, a cursory look suggests 
that SFGs appear to form stars with an S-curve age dis- 
tribution. In that case, average present-day SFGs grew 
rapidly at high redshift (z > 2) before biding their time 
at intermediate redshifts - while sSFRs were still 3—10 
times higher than today. Finally, they surge again, form- 
ing stars rapidly over the last 1 — 2 GyrQ Such present- 
day SFGs would have occupied highly biased position in 
the SFR main sequence, possibly leaving their mark in 
environmental or structural trends. 

However, such a naive interpretation of the data is un- 
warranted in light of uncertainty in stellar population 
ages. This uncertainty drives the peak of SFR to spread 
out, and the logarithmic nature of variations in stellar 
populations naturally gives rise to that S-curve in cumu- 
lative star formation. Specifically, intermediate to old 
populations shift mass from a high(er) redshift SFR peak 
to low redshift, while a tail of star formation extends to 
the beginning of the universe (see figures in Appendix [B] 
also). Whatever the SED or CMD quality, the nature 
of this shape change should be a generic part of prop- 
erly regularized SED-based analysi£3. The impact on 
"true" MSI predictions (blue) is a shift to SFHs that are 
well matched to the VESPA lBCOrH 1 ! SED analysis of SDSS 
(magenta) at all masses shown. 

4.4. Star Formation Trends and Local Dwarfs 

For the I Will analysis of dlrr, on the other hand, Fig- 
ure [7] shows significant differences from extrapolation of 
MSI below the mass completeness range of SFR surveys. 
Our simp le mo del of age errors and those errors esti- 
mated in I Will are unable to account for these differ- 
ences. In this regime, MSI predicts excessive downsizing 
such that dwar fs form most of their mass below z = 1, 
while the I Will dlrr galaxies form most of their mass at 
z > 2. 

11 The cumulative star formation fraction is more easily com- 
pared to observations than M«/M»o- However, the two are similar 
and would be identical in the case of fixed mass return fraction 
(i.e., constant R(t)). 

These fe atures were noted in past analysis lIToieiro et al.1 2009; 
IPanter et aLll2007l1 , but attributed to TP-AGB stars because IM05I 
shifts some of the low redshift star formation to earlier epochs 
(although the S-curve shape is still visible). 

3 Although the proper error model may vary depending on fea- 
tures in SSP spectral templates. 

14 The match to BC03 as opposed to M05 should be interpreted 
cautiously. The SFHs reflect some combination of the real inversion 
of flux from the SED, and the way noise is interpreted in terms of 
the SSPs basis. A log-age error is expected, but it may not capture 
all important artifacts in the flux inversion. 
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Figure 6. Cumulative star formation from the SDSS SEP modeling (poin ts) and MSI (shaded regions) in four panels for different mass 
samples. VESPA's SED analysis of SDSS is shown for the M05 (triangles), and BC03 (squares) population synthesis models. MSI (blue) was 
initialized with the same average mass and redsh ift a s tha t of the SDSS samples and then smoothed to mimic age errors (magenta solid). 
The shading encompasses variation between Kll and O10 SFR main sequence observations. Striping indicates the reliability of the results 
as described in Figure [3] 



Figure [5] demonstrate s the e xtent of this shift in the matched by MSI trends. The figure shows SFHs in 
context of the VESPA IBC03I SED analysis that was the range 1O 11 - 1L5 M , 10 10 - 10 - 5 M©, 10 9 - 9 5 M Q , and 
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Figure 7. Same as Figure \E\ but black points arc Wll CMD- 
based results for dwarf irregulars (circles). Error bars show sta- 
tistical uncertainty on the mean, which is similar to systematic 
variation between different stellar evolution models. MSI trends 
are extrapolated over this mass range. 



1O 8 " 8 5 M (from lightest to darkest), plotted over the 
lower mass IW ill dlrrs (most with M* < 10 8 M Q ) that are 
scaled to the linearly interpolated dwarf mass at (z b s ) of 
the given sample. Archaeological downsizing is apparent 
in the p urely star-forming sample (as previously note d 
by, e.g., iCid Femandes et all 120051: lAsari et all 120071 ) . 
Nevertheless, the dlrr galaxies seem to reverse this 
trend: they assemble more like galaxies of ~ lO n M0, 
and deviate str ong ly from their ~ 1O 8 M siblings ana- 
lyzed with VESPJo an d trends in the SFR m ain sequence. 

This is in line with a central conclusion of lWeisz et al.1 
(|2011al) . that the formation of dwarfs in the local vol- 
ume agrees with the cosmic SFH (which is dominated by 
massive SFGs, and may be overestimated at z > z orr ), 
but here we have highlighted that these lowest mass 
dwarf SFHs seem quite inconsistent with the cosmologi- 
cal downsizing trends found in SFR observations and the 
SED-based fossil record. A discussion of this shift, and 
its observational context follows in 35.21 

5. DISCUSSION 

In section [3] we showed that MSI indicates late for- 
mation of typical SFGs (z < 2 and tvo S; 10 Gyr), then 
compared the corresponding mass growth with the fossil 
record in section |U In this section we place these MSI 
results in a broader physical and observational context. 

15 To ensure that the downsizing trend we present at low masses 
was not purely the result of noisy SEDs, we also limited our sample 
to only galaxies with high quality spectra with average S/N> 20 in 
the regions that are fit ((S/N) = 27), and found identical results. 



Figur e 8. Cumulative star formation for SDSS SFGs from VESPA 
IBC03 SED modeling in mass bins at lO 8-85 , 10 9 " 9 ' 5 , 10 10 ~ 10 ' 5 , 



1 gio.5— li ^ anc [ io 11— 11-5 ) colore d (squ ares) from lightest to dark- 
est respectively, compared to the IW11I dlrr sample (black circles) . 
Each mass bin has been rescaled to the linearly interpolated cu- 
mulative star formation fraction of the dlrr sample at their epoch 
of observation. Star formation moves to lower redshift as mass de- 
creases (downsizing), but the dlrrs form more like the most massive 
sample. 



5.1. A Simulation Context 

Main sequence integration has no connection to galaxy 
dynamics, so it is interesting to note that its implications 
align with the trends seen in increasingly successful cos- 
mological simulations. In particular, the late-formation 
of SFGs suggested by MSI may resolve the (potentially 
orthogonal) struggle to reproduce the structure of late- 
type SFGs in simulations. 

Current galaxy formation simulations typically con- 
vert too much low angular momentum gas into 
massive bulges. The observed median bulge-to- 
total (B/T) stellar mass ratio is 15% in bright (> 
10 10 M Q ) local late-ty p e gal axies (jWeinzirl et al.1 120091: 
IGadotti fe Kauffmannl [2001 . and "bulgeless"E3 disk 
galaxies are abundant in the local universe (espe- 
cially dominating at < 10 10 ' 5 M Q iKormendv et alj|2010t 
iFisher fc Drory||201lD . Such small bulges, embedded in 
star-forming disks, are rarely reproduced in the current 
generation of simulations of M*o > 10 10 M Q galaxies. 

Importantly, whatever the physical or numerical 
mecha nisms that drive s i mulat e d bulge formation 
(e.g.. IPiontek fc Stemmed [2001 IBrook et all [20TT1: 
iGovernato et al.l I2010D . state-of-the-art simulations of 
SFGs consistently produce their excessively massive 

16 A bulgeless g alaxy is defined as having Scrsic index n < 2 
(see, [Graham 20111, for further discussion). 
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bulges at 2 > 2 - earlier than MSI indicates ~ 
15% of stellar mass should have assembled. S pecifi - 
cally. iScannapieco et all (|2011l ) and lStinson et al.l (|2010D 
find that their kinematic bulges comprise the major- 
ity of each one of their 17 galaxy's stellar masses; the 
Scannapieco et all (20111 ) bulges form at z > 2, and the 



Stinson et al.l (|2010t Figure 14) non-disk stars a lso form 



early. The higher resolution iBrooks et al.l ([2011D sample 
of 6 massive galaxies with B/T > 0.26 ((B/T) = 0.44) 
also form at least the majority of their mas s at z > 1 (A 
Brooks private communication; also see, IBrooks et al 

mand discussion of sSFRs therein). IGuedes et al 
simulated a Milky Way-like disk galaxy with a 
(photometric) B/T = 0.26, again, with almost all of the 
bulge mass forming at z > 2. Suppressed star formation 
at z > 2 is responsible for the reduced bul ge fraction in 
a para meter stu dy of a very massiv e disk in lAgertz et al.l 
(|20ll . Finally, IBrook et al] (|20TTI ) forms a B/T = 0.21 
disk with a SFH that closely matches our fiducial SFH 
expectation for a 10 10 M Q system (Figure |3|). Delayed 
formation seems not only achievable in simulations, but 
also may be important for reproducing typical disk struc- 
tures. 

Intriguingly, two sm all dwarf galaxy simulations by 
iGovernato et al.l (|2010l ) manage to form early without 
forming massive bulges (one is fit by an exponential pro- 
file and one has bulge-to-disk ratio of 0.08). Such simu- 
lations may be important for understanding suppressed 
star formation at low masses and high redshift for all 
galaxies. They may also l ead to a resolution of the unan- 
ticipated early growth of IW11I dwarfs (but see the next 
section) . 



5.2. Resolved and High S/N Observations 

The a brupt shift in the CMD-based age distribution 
of rWTll dlrrs compared to SED- and MSI-based SFHs, 
could instead highlight problems with MSI-, SED-, or 
C MD-based SFHs. T he handful of massive galaxies 
in iWeisz et al.l (|2011bl) form even more of their mass 
late, with the most massive galaxies in their sample 
(M* > 10 7 - 75 M o ) reported to have formed 64 ± 0.6% of 
their stars at z > 2. That is slightly more z > 2 star for- 
mation than in the lower mass dlrrs (52+g%) and makes 
the nuances of dwarf galaxy formation a less persuasive 
explanation for the transition to early-forming dwarfs. 

Another explanation is simply that the sample is small 
and biased by t he local group overdensity. However, 
those same IW 111 dwarfs fit nicely into sSFR trends when 
measur ed as part of the larger G ALEX-HUGS 11 sam- 
ple of iLee et al.1 ([2009B)(sSFRkii [10 7 6 M q ,z = 0] « 
sSFRLeeuv^O^Gyr" 1 ). 

Unfortunately, other analyses of resolved populations 
in SFG disks are limited in number and ability to dis- 
cern age differences, with too-large error bars on mass 
formed in a given bin for statis tically significant infer- 
ences fe.g.. lWilliams et al.ll20rl . Where any constraint 
on the z = 2 divide is reported, results are also mixed: 
NGC 300 (jOogarten et alJl2010l . 2.4 x 1 10 M^) formed 
70% of its stars at z > 2, while M33 dWilliams et al ' 



2^0 llBarker"et^l2007l) . and NGC 2967 (IWilliams et al 
20101) may have formed m ost of their mass at z < 2. 
iMacArthur et all (|2009D and iSanchezT Slazqucz et al.l 



(|2011l) recently performed full spectrum fitsQ to high 
S/N (> 50) data in the central regions of nearby disk 
galaxies. While these measurements target the inner 
disks and bulges of galaxies, the results endorse strikingly 
old populations for the bulk of stellar mass in galaxy 
disks. Of their collective 10 galaxies (M* ~ 10 Mq) 
with measurements to the disk scale- length (r^), six 
(10239, N7495, N7490, N0173, N1358, and N1365) ap- 
pear to form most, in several cases more than 90%, of 
their stars at before z = 2. The blue regions in Fig- 
ure [6] shows that galaxies of this mass would be expected 
to form 10% of their mass over the last 5 Gyr. If such 
galaxies form 90% of their entire stellar mass at z > 2, 
they would have persisted as massive quenched disks for 
around 5 Gyr. 

If born out as representative, these results may demand 
either an important recalibration of SPS models, or a 
shift in our understanding of where present-day SFGs sit 
in the SFR main sequence (something not apparent in 
studies of environment, or the structural main sequence, 
see £13. 1 ,2|) . More generally, such extreme examples high- 
light the value of cross checks between the SFR main 
sequence and fossil reco rd analysis, e s pecially of the sort 
recently undertaken bv lWuvts et al.l (|2011bl ). Clearly, a 
larger number of high-S/N observations and robust ex- 
aminations of disk ages will be crucial for a better under- 
standing of SFG growth and of the galaxy fossil record. 

6. SUMMARY AND CONCLUSIONS 

We have derived the average growth of stellar mass 
in present-day SFGs by simple integration of the consen- 
sus SFR main sequence under the assumption that SFGs 
were always SFGs in the past (an approach that we term 
MSI). Our results regarding average stellar mass growth 
can be summarized as follows: 

• SFGs form late such that only ~ 15% of the stellar 
mass in the progenitors of 5 x 10 9 M Q < M*o < 
1O 11 M0 galaxies was in place before z = 1 — 2 
(depending on mass and downsizing). 

• It follows that massive SFGs at z > 2 are not ex- 
pected to be progenitors of typical SFGs today. 

• The effect of mergers on MSI results was found to 
be negligible based on semicmpirical merger trees 
generated in the MSI framework. 

• Similarly, there is no clear evidence that the way 
SFGs occupy the distribution of SFRs in the 
SFR— M* relation significantly alter MSI results. 

• Accurate analytic approximations to average SFHs 
derived from MSI are presented in Appendix [A"l 

Moreover, the delayed formation of SFGs implied by 
our analysis (but f irst noted in the stage-T models of 
iNoeske et a"ni2007al) was found to be consistent with ob- 
servations of the evolution of stellar mass in the uni- 
vers e, is also implied by halo-based semiempirical meth- 
ods (jConrov fc Wechslerl 12001 1 Zheng et alJl2007aD . and 

17 At a minimum full spectrum fitting with multiple SSPs is 
critical to age determinations of SFG because of their extended 
SFHs. Single SSPs fits or Lick indices by themselves are insufficient 
as detailed in UMacArthur et al.|[2009| ) 
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fits well in the context of increasingly successful cosmo- 
logical simulations. 

When comparing MSI-bascd SFHs to those inferred 
from the fossil record, we found that, 

• Expected age uncertainties in SED-based analysis 
cause a characteristic shape distortion in derived 
SFHs. 

• SED analysis of SDSS SFGs could be reconciled 
with MSI after accounting for these expected age 
uncertainties. 

• Local dwarf galaxies with (M*) < 10 s , on the other 
hand, were found to depart from MSI downsiz- 
ing extrapolations, with SFHs resembling those of 
10 11 M Q galaxies. 

We stressed that the last point may reflect systematic er- 
rors in population synthesis models rather than a phys- 
ical transition. For a handful of cases, CMD and high 
S/N SED analysis of more massive galaxies suggests more 
early star formation than inferred from MSI. A larger 
sample of SFGs is greatly needed and it is clear from our 



SDSS comparison that a full understanding age uncer- 
tainties is crucial for interpretation of that data. It will 
be interesting to see how improved SED modeling and 
larger data sets can be accommodated in the SFR main 
sequence. 
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APPENDIX 
ANALYTIC SFHS 

With a few approximations, accurate analytic formulae for average SFH and stellar mass growth based on SFR 
observations can be derived. Since these may b e useful for tying chemical evolution or SPS models to SFR main 
sequence measurements (e.g. JWuvts et al.ll2011bO or for simulation comparisons, they are provided below. 

First, in the concordance cosmology, the expansion factor is linear as a function of lookback time to < 5% accuracy 
until iib « 12.5 Gyr or z w 4, 

(1 + z) « (1 - A ii b )-\ (Al) 



where A a = (1 — a(£ii n ))Aiin and we take tn 



12 Gyr so that A a = 0.064 Gyr 1 in our cosmology. Rewriting Eq. [3] in 

3i(t)/<I>(t) being the fraction of the current 



terms of the stellar mass fraction formed = , and with R(t) 
SFR returned to the ISM by stellar mass loss, then, for a power law fit to the SFR main sequence (Eq. [T]), 



where /o 



A/.o 



hit) 

Taking R(t) 
/*(*) = 



[l-R(t)]foMt) 1+p (l~\J) 
0.45 and integrating for /*, 

13 fo [1 - R(t)] 



(a - 1) A a 



[(l + z)- 1 -!] 



(A2) 



(A3) 



Then, plugging Eq. IA3l into Eq. IA2[ analytic SFHs, $(i) = M*o/*(i)/ [1 — can be synthesized from a power law 

jr. 



$(t) = M* f (l + z) c 



(a - 1) A n 



1 + 



(A4) 



a and j3 parameters can be found in Table [T] (under the "pow" form) and are commonly reported with SFR sequence 
observations. Any broken power law in a or (j is a simple extension (e.g., for a flattening of a at z ss 2 or a flattening 
of (3 in dwarfs). 

The dashed gray line in Figure |9] shows analytic SFHs and /*(t) for constant R(t) = 0.45. Results are in good 
agreement with the full MSI procedure, but late mass growth is overestimated in massive galaxies because mass los s 
from old stars is not properly accounted for when assuming a constant return fraction (see iLeitner fc Kravtsovll201ll) . 

As a first order correction for massive galaxies, the return fraction can be decomposed into a constant term from 
young stars, Rq, and a term related to accumulated old stellar mass. The relative fraction of gas returned by old stars 
is roughly proportional to a galaxy's age (if or m — hb) (i-e., amount of stellar mass if star formation were constant since 
iform) calibrated to a reference return fraction, i?n, in a galaxy of lO n M0 with iformii- Then 
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Figure 9. SFHs (left) and mass g rowth (right) from MSI, in galaxies of M, = 10 s M , 10 9 M Q , 10 10 M Q and 10 11 M Q . Calculations 
are shown for the fiducial SFR fit l IKHI ) using t he fu ll MSI approach (blue lines), and analytic approximations with 5R = 0.45 (gray 
long-dashed), and 5ft (i) (black dotted; given by Eg. IA5II . Vertical dotted lines are z eTT from Figure [2] 



assigning /*form(*form) =0.1: 



^a^form 1 



(a-l)A a (l-/;£ m ) 
P fo [1 - R(t)} 



(A6) 



We calculate tf OTm values using R(t) = 0. 45, and plug those values into Eq, IA51 with i?o = 0.45 and R\i = 0.30. 
Inserting the full return rate from Eq. IA5l into Eg. I A3I results in the <!>(£) and /*(£) values plotted (dotted lines) over 
the full MSI models (solid lines) in Figure El The /* and SFH approximations are accurate to a couple of percent of 
the final stellar mass (i.e., 20% error when the galaxy is 10% of M*o)- The same correction are equally accurate for 
the lOlOl parameterization. 

NOTES ON AGE RESOLUTION 

In q4.3l we show that SFHs based on the fossil record can be misleading, particularly about early galaxy growth. The 
root of this cl aim is the consensus view that the max i mum reliable precisi on of age determinations in SED-analysis is 
low (> 0.5 dex lOcvirk et alJfeOOfl iWalcher et aT1l201 It [Cid Fernandesll2007l ). Our solution is a similarly uncontroversial 
smoothing of the SFH to mimic age errors. The result, however, is a notable change in the shape of SFHs that brings 
VESPA results into line with cosmological SFR trends (i.e., the SFR main sequence). I n light of the import ance of noise 
modeling to this result, a discussion of age errors is prudent. For guidance, we rely on lOcvirk et al.l (|2006| ) who carried 
out a thorough simulat ion campaign to quantify the noise properties of SSPs and full spectrum fitting in PEGASE-HR 
(|Le Borgne et al.l 120041 ) with implications for interpretation of any optical SED inversion^. This section provides a 
short summary of the issue, and a discussion of how it pertains to VESPA results. 

Modeling the intrinsic SEDs of galaxies involves the conflicting goals of (1) finding a stellar population (metallicity 
and extinction) that generates a good fit to SED, and (2) not letting that population be shaped by noise in SEDs. 
Even fitting in the absence of noise is not clear cut because stellar evolut i on models do not match a l l of the observed 
featu r es in galaxy spectra (see e.g.. iToieiro et al.ll2007t IPanter et al"1l2007t ICid Fernandes et aT1l2005t [Conrov fc Gunnl 
120101: IWalcher et all 1201 ID and the calibration of i mportant pieces of stellar evolut ion (e.g., TP-AGB stars at long 
wave l engths: lMaraston||2005 ; lConrov k, G mm 2010; Maraston 2011; Kriek ct alj|2010l : stellar abundance: IPanter et al.l 



120071 iMacArthur et al.ll2009t IMF, etc.) are ongoing. lOcvirk et al.l (|2006f ) found that noise calibrated into the model 
PEGASE-HR spectra themselves place a fundamental floor of ~ 0.4 dex on age recovery for even the highest S/N 
spectra. They also found that noise in SEDs typical of SDSS spectra (S/N« 20 per ~ 3 A pixel in the unmasked VESPA 
range) renders differences between simple SSPs unrecoverable within ~ 0.5 dex. Any smaller-scale time information 

18 | Ocvirk et al. | (2003) conclusions were reached over the slightly 3SQ0Q _ 7gooA (including masked emission Unes) . 

smaller spectral range, from 4000 — 6800A compared to VESPA's 
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Figure 10. SFH (left) and cumulative star formation fraction (right) for SDSS SFGs of 1O 1O - 5_11 M0 from VESPA analyzed with 
Bruzual & Chariot (2003) (solid black), compared to MSI (dashed) with varied age resolution (see legend). MSI starts with the same 
average mass and look back time as the SDSS sample. The first 7 VESPA age-bins were combined to remove artifacts. 



derived from SEDs must come from noise, improved modeling, or the use of a larger spectral range. Moreover, these 
resolutions are further broade ned when considering non-SSPs with uncertain mctallicities, extinction and physical 
effects (e.g.. rWuvts et aLll2009D . Regularization techniques can be used to assign resolution to solutions so that they 
are consistent with these noise limitations and primarily responsive to signal. 

VESPA attempts to regularize the stellar population solution using the general procedure of lOcvirk et al.l (|2006l 
STECMAP), but adaptively refined logarithmic age bins for basis vectors. Given a solution basis of age bins, metal- 
licities and extinction, VESPA estimates the contribution of noise to the coefficients (parameters) that multiply the 
solution basis. The number of p arameters used to d escribe the solution should be < re, the number of coefficients 
that are not noise dominated (see lToieiro et al1l2007l §2.2.2 ). Then, to select an appropriate age basis, the bins that 
contribute the most to the total flux in the model are iteratively refined until they reach the maximum A agc « 0.2 dex 
resolution, re is recalculated at each iteration. Finally, the chosen solution is the most refined solution that has fewer 
non-zero flux age bins (and o ther model parame ters) than re. 

Contrary to the results of (jOcvirk et al.ll2006l ) that A age < 0.4 dex is not achievable, VESPA regularly refines to its 
highest resolution A age « 0.2 dex. A full exploration of this discrepancy is well beyond the scope of this paper, but one 
plausible explanation is that, since bins containing no flux do not count toward the parameter limit, re, VESPA favors 
higher frequency solution variations than can be robustly extracted (R. Tojeiro, private communications). In contrast, 
STECMAP preferentially weights smooth basis vectors for their solutions using well tested algorithms to elim inate 
artifacts. Although it is not clear that the smoothed STECMAP basis is ideal (as noted bv lOcvirk et al.lf2006f ). the 
VESPA basis appears to be to aggressive. Ill conditioning of the inversion means there is no way to know exactly how 
noise-modified light gets interpreted into a stellar age spread. 

Artifacts are to be expected but, supported by covariance primarily between adjacent bins in VESPA tests, we assume 
that VESPA will mimic the STECMAP regularization. That regularization favors small gradients in log(£ age ), finds 
little bias in (log(i age )}, and almost constant logarithmic resolution ±A agc /2 about any median log-age value. In this 
case a log-age Gaussian convolution with constant full width ha lf max (A ayR ) is the obvious choice for a noise-model. 
The 1 dex resolution, as noted in §4.21 is motivated both by the (j Tojeiro et al.l [2007] ) covariance between adjacent bins 
and the typical mass- weighted size of the bins. The convolution preserves (log(i age )) and biases (i age ), depending on 
how the iago distribution is limited by cosmology. 

Figure [TO] demonstrates that the A age size for massive SFGs, does not qualitatively alter our results. The figure 
shows the SFH and cumulative star formation recovered by VESPA for galaxies of 1O 1O_1O 5 M0 using iBruzual &: Chariot] 
(|2003l) analysis binned and averaged at the highest age resolution (black solid line) , compared to MSI models (dashed) 
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smoothed in log-age with A ag0 = 0,0.5,0.75 and 1 dex. The first seven VESPA bins were combined to hide (dramatic) 
artifacts in these bins The S-curve noted in £14.31 grows (right) with increasing A agc and, correspondingly, the SFH 
peak flattens (left). 
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